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Engine  Research  Center 
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ABST^T 


High  injection  pressure  impinging  spray 
experiments  and  modeling  were  performed 
under  simulated  diesel  engine  conditions 
(pressure  and  density)  at  ambient  temperature. 

A  spray  impinged  normal  to  a  small  crown  in 
the  bowl  of  a  simulated  piston.  High  speed 
photography  was  used  in  the  constant  volume 
bomb  to  examine  the  effect  of  impingement  on 
fuel  mixing.  The  spray  model  which  includes 
drop  breakup,  coalescence,  impingement,  and 
vaporization  effects  was  used  to  predict  fuel 
mixing  in  the  bomb.  The  spray  distributions 
predicted  by  the  model  are  compared  to  the 
photographs  obtained  in  the  bomb.  * 

SPRAY  MPWGEMENT  BACKGROUND  er) 

Interest  in  direct  injection  stratified  1 
charge  engines  and  direct  injection  diesel 
engines  has  included  strong  interest  in 
developing  accurate  models  for  liquid  fuel 
sprays.  Among  the  more  widely  used  models 
is  KIVA,  a  code  developed  at  Los  Alamos 
National  Laboratory  f  t  ]*.  Recent 
developments  have  focussed  on  modeling  basic 
gas  phase  fluid,  mechanics  and  using  liquid 
spray  submodels  to  describe  the  behavior  of 
the  liquid  spray.  The  detailed  spray  submodels 
have  primarily  focussed  on  improved  modeling 
of  drop  breakup  [  2,3  ],  drop  ooatoscence  [  4  ], 
and  improved  toteoriptions  of  the  initial 
conditions  and  bounctory^eondftfons  applied  to 
the  problem  (  4  ].  The  results  of  this  work  have 
been  evaluated  by  comparison  with  limited 


an? 


amounts  of  data  for  sprays  in  regions  far  from 
toe  injector  (>200d  where  d  is  the  injector 
orifice  diameter).  These  results  have 
emphasized  the  significance  of  the 
combination  of  drop  breakup  and  coalescence, 
and  the  relative  insensitivity  of  far  field  drop 
size,  velocity,  and  spread  angle  results  to  the 
injection  conditions.  In  addition  to  studies  of 
liquid  sprays  in  a  semi-infinite  gas  field, 
model  results  have  been  extended  to  finite  gas 
fields  where  the  spray  impinges  on  a  solid 
surface  at  some  location  relatively  far  from 
the  injector  tip  (  >200d).  A  recent  paper  by 
Naber  and  Reitz  [  5  J  has  addressed  this  issue 
by  examining  the  problem  of  a  diesel  spray 
impinging  on  the  bowl  of  a  piston  in  a  moderate 
swirl  engine.  Comparison  of  their  results  with 
photographs  from  an  operating  engine  indicate 
good  qualitative  agreement  between  their 
model  predictions  for  overall  far-field  spread 
and  apparent  surface  wetting  and  the  data 
acquired  from  toe  photographs. 

In  some  combustion  situations  of 
interest,  the  geometry  of  the  combustion 
chamber  may  limb  the  extent  of  the  spray  to 
regions  of  significantly  smaller  i/d  values 
than  those  for  which  model  results  have  been 
established.  In  these  near  field  regions,  the 
issues  of  spray  behavior  near  the  nozzle  tip 
and  the  effects  of  high  velocity  liquid  sprays 
impinging  on  solid  surfaces  may  be  aigniftoant. 

Proposals  for  relative^  novel  cylinder 
geometries,  In  particular  novel  piston 

e  ia  «edM«^yhaa 
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,  mod  ot  the  end  of  this  paper. 
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QMffto tries,  have  recently  been  made  by 
Kroeger  [  6  ]  and  Kato  and  Onishi  [  7  ].  These 
obambe*  detlgns  generally  include  a  piston  t 
bbh  a  boM  arid  a  smaR  projection  in  the 
mkkU*  it  the  bowl.  These  piston  top 
geometries  appear  'w  have  been  developed  for 
somesrhft  different  applications:  Kroeger's 
(Caterpillar)  for  direct  injection  diesels  using 
neat  methanol,  and  Kato  and  Onishi's  (NICE)  for 
direct  injection  stratified  charge  (DISC) 
engines.  The  basic  geometry  exhibits  a  raised 
surface  at  or  near  toe  center  of  the  piston 
bowl,  onto  which  some  or  all  of  the  fuel  spray 
impacts.  In  the  case  of  the  Caterpillar  engine, 
nine  other  spray  lobes  are  also  injected  in  a 
normal  pattern  radially  and  slightly  downward 
from  toe  injector.  In  the  case  of  too  MCE 
engine  alt  of  the  fUetis  sprayed  onto  toe  raised 
piston  impingement  area  from  a  single  orifice. 
In  either  case  the  liquid  which  impacts  the 
raised  surface  is  axpebtsd  to  provide  a 
centrally  located  cloud  of  fuel  drops  and  fuel 
vapor  which  will  enhance  the  subsequent 
combustion  event.  In  both  cases,  toe  intent  of 
the  piston  projection  is  id  enhance  the  fust 
injection  process  by  directing  tome  or  ell  of 
the  fuel  to  be  injected  onto  the  projection. 

This  projection  is  typicaHy  15-20  nozzle 
orifice  diameters  downstream  of  the  nozzle 


tip,  which  puts  it  in  a  region  of  toe  very  near 
spray  field  where  liquid  velocities  are  quite 
high.  Improved  performance  or  continued 
performance  under  atypical  operating 
conditions  were  demonstrated  experimentally 
in  both  of  toeee  engine*.  In  general,  this  is  a 
very  different  situation  from  the  conditions  of 


previous  work  in  spray  modeling,  yet  it 
presents  an  area  of  particular  Interest  in 
terms  of  enslyitrtg  toe  petfewnanos  of 
aigntoeantly  different  combustion  and  fuel 
injection  geometriSe.  These  new  geometries 
may  present  a  parifeutor  challenge  to  odetorg 
codes,  aw  toev  toctode  efrrdcsi1  ath&m  hot 

^ wl  fBBqj  SUCH 

wPWwW  WP  wIMVlMWIft  Vllr  a  WmUQWm' 

: MUMuaediMtoiiiatfliriifbi  **-■  -  —  - •—  -» 

lw  wf  eO®  I^WSr  wQW  wf 


This  paper  presents  the  application  of 
one  currently  available  spray  code  (KIVA)  to 
toe  situation  of  toe  piston  geometry  displayed 
by  the  engines  presented  by  Kroeger  {  6  ]  and  by 
Kato  and  Onishi  ( 7  ].  A  set  of  experiments 
has  been  performed  for  the  raised  piston 
center  section  geometry  using  a  high  pressure 
diesel  fuel  injector  in  a  high  pressure  unheated 
bomb.  High  speed  movies  provided  information 
on  the  spread  and  location  of  the  spray  under 
different  operating  conditions.  The  results 
from  the  movies  are  compared  to  toe  results  of 
the  computer  model  in  its  original  state,  and  in 
a  modified  state. 


NUMERICAL  MOOEL 


A  three  dimensional  finite  difference 
numerical  model  (KIVA)  was  used  to  model  the 
spray  in  toe  constant  volume  bomb.  The  model 
aolvaa  toe  gaa  phase  using  a  time  explicit 
euierian  finite  difference  scheme  for  the 
averaged  Navier-Stokes  equations  with  a  k-  e 
turbulence  model.  For  this  application,  the 
geometry  of  the  spray  was  considered 
axisymmetric,  so  a  two-dimensional  version  of 
toe  code  was  used.  The  spray  is  modeled  using 
a  stochastic  lagrengian  scheme  where  parcels 
represent  a  number  of  drops  with  common 
properties.  The  parcels  interact  with  toe  gas 
through  terms  in  the  gas  equations  exchanging 
mass,  momentum,  and  energy,  and  generating 
turbulence.  Drop  coalescence,  breakup  and 
impingement  are  included  in  the  stochastic 
parcel  model.  More  details  of  the  numerical 
model  are  given  in  toe  references  (1 ,2,5,8]. 

The  breakup  model  used  in  this  study  is 
described  in  detail  by  Reitz  and  Diwakar  [  8  ]. 
The  initial  parcel  diameter  at  the  time  of 
injection  it  equal  to  the  nozzle  diameter.  This 
differs  from  tiieir  earlier  work  [  2  ]  and  that  of 
Ameden  et.al.  [  1  ]  where  the  initial  parcel 
diameters  were  determined  from  distributions. 
Ametierfs  (  1  ]  probability  distribution 
rtomitfizatf  to  (toft  total  weight  is 


KD 


'  e% 


%r  4 
**» 


•  (1) 


Where  rM  is  toe  Sauter  mean 

’  tnUftthn  ikA 

I  the 

WM  m  ww  BfMHS  W 


881316 


3 


coalescence  determines  the  drop  distribution 
downstream  of  the  injector.  Two  modes  of  drop 
breakup  bag  and  stripping  are  modeled.  Bag 
breakup  occurs  for 


and  stripping  breakup  tor 


(2> 
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(3a, 3b) 


The  impingement  of  drops  on  solid  boundaries 
is  modeled  using  an  analogy  with  potential 
flow  jet  impingement  as  described  by  Naber 
and  Reitz  [  5  ].  The  drop  leaves  tangent  to  toe 
surface  in  a  direction  determined  from  a 
probability  distribution  function.  Hie 
probability  distribution  function  is  derived 
from  potential  flow  jet  conservation  of  mass 
and  momentum  and  an  assumed  momentum 
distribution.  The  angle  ty  is  the  angle  of  the 
tangential  velocity  of  toe  parcel  after 
collision  with  the  wall  in  the  plane  of  the 
surface,  relative  to  the  impinging  tangential 
velocity  vector.  The  parcel’s  new  velocity 
vector  direction  is  determined  from  the 
equation 


y=-£ln{l-x(l-e1} 

P  (4) 


where  X  is  a  uniform  random  number  (0,1),  and 
p  is  determined  from  the  function 


m(a)  * 


e^+1 
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Here  a  is  the  in* tngement  angle  measured 
wmoi  neonm  rnepereoi  oosanot 


due  to  the  wan 


conditions  and  boundary  conditions.  For  toe 
gas  flow,  the  initial  condition  is  a  uniform 
field  of  a  known  pressure  and  density.  The 
initial  gas  velocity  is  zero  and  the  temperature 
is  a  constant  throughout  the  problem,  as  the 
bomb  is  operated  at  room  temperature  and 
vaporization  is  neglected  in  toe  model.  For  toe 
spray,  an  initial  drop  size  or  distribution, 

Initial  velocity,  and  spread  angle  must  be 
specified.  As  toe  initial  drop  size  was  varied 
during  toe  course  of  the  work,  discussion  of 
the  actual  values  used  is  included  in  the 
discussion  of  toe  model  results.  Exponentially 
decaying  spray  velocity  profiles  were 
estimated  using  the  known  peak  (initial)  fuel 
pressures,  injector  closing  pressures, 
injection  durations,  fuel  delivery,  and  orifice 
discharge  coefficients  provided  by  the  injector 
manufacturer  [9].  The  initial  spray  angle 
was  chosen  from  the  photographs  of  the  spray. 
The  grid  used  in  the  computations  is  shown  in 
Fig.  1.  A  single  azimuthal  sector  of  0.5 
degrees  was  used  for  the  axisymmetric 
geometry  with  17  radial  and  13  axial  grid 
points.  A  typical  run  with  this  model  took 
about  30  minutes  of  computer  time  on  a  Cray 
)WP. 

For  each  set  of  conditions  modeled,  toe 
output  consists  of  a  series  of  plots  at  various 
times,  providing  the  current  location  of  the 
spray  parcels,  liquid  and  vapor  concentration 
contours,  and  temperature  contours.  In  this 
paper  only  the  spray  parcel  location  plots  will 
be  displayed. 


Ffg.l.  Axisymmetric  computational  grid. 
Domain  used  far  calculation  of  impinging 
•pray  to  constant  volume  bomb. 
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Front  View 


Fuef  Injector  Mount 


Section  A-A 


Fig.  2.  Front  and  section  views  of  pressurized  bomb  with  optical 
access.  End  caps  orientated  to  view  top  portion  of  chamber. 


CONSTANT  VOLUME  BOMB  EXPB*M0VTS 


A  constant  volume  bomb  was  constructed 
to  provide  a  chamber  capable  of  providing  the 
high  gas  pressures  and  densities  common  in 
high  compression  ratio  engines.  A  sketch  of 
the  constant  volume  bomb  is  given  in  Fig.  2. 

The  bomb  was  constructed  of  203  mm  (8  in.) 
diameter  pipe,  203  mm  (8  in.)  long.  At  each 
end  of  the  pipe,  a  flange  was  welded  on  and  a 
cap  fitted  to  the  end.  In  the  caps  on  each  end 
were  mounted  101  mm  (4  in.)  quartz  windows. 
The  window  arrangement  was  designed  to 
allow  the  and  caps  to  be  remounted  bt  a  variety 
of  positions  to  allow  a  view  of  various 
portions  of  the  spray.  At  the  osntsr,  axially,  of 


the  chamber,  a  mounting  ftonoe  was  welded 
into  the  eWf  we#  for  mounting  the  fuel 
injector.  T1iefpoaMOn  M  injector  and  tip 
warn  acjjustoOfo  fo  jomi  exfont  to  ensure  the 
shape  of  the  vofcpto  pfoyett  Kttiw  rote  in  the 
suhMouMt  aonv  rm^qa  Thi  chimbor  hit 

wwirvwi|iivin  aynam  itnrmrtf*  r  rw  unminrvi  inw 
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single  hole  0.406  mm  diameter  nozzle  and  a 
maximum  peak  fuel  pressure  greater  than  138 
MPa.  This  injection  system  exhibited  a  rapid 
needle  lift  to  introduce  the  fuel  at  a  pressure 
near  the  peak  fuel  pressure.  Peak  pressure  was 
reached  almost  instantaneously,  and  dropped 
steadily  until  a  preiced ed  spring  ended 
injection.  The  quantity  of  fuel  injected  is 
determined  by  the  injector  geometry  and  the 
fuel  rail  pressure.  Fig.  3  shows  a 
representative  injection  pressure  profile  taken 
from  [  9  ]  for  this  type  of  injector. 

The  pressure  and  density  in  the  bomb 
were  chosen  to  match  that  of  a  test  engine  at 
15  degrees  btdc  (P0-32.0  atm,  p9-14.9  Kg/m3  at 
T-760K)  by  using  a  mixture  of  0.30  N2  and  0.70 
He  by  volume.  Table  1  lists  injector  conditions 
for  the  three  events  photographed.  The  fuel 
ueed  in  all  of  tne  experiments  was  standard 
type  I  referee  fuel. 

*Phe  time  varying  extant  of  the  fuel  spray 
was  measured  in  the  bomb  using  a  high  speed 
movie  camera  running  it  5000  frames  per 
•8bbfl&  A  oopper-vapor  laser  providing  pulses 
approximately  to  na  long,  at  fppitllion  rates 
up  to  10  kHz  was  ueed  as  a  light  source.  For 

stswfwettizeti  to  Mill  ewe  te&m  pule#  ppr 
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Fig.  3.  Typical  injector  pressure  profile 
of  BKM  injector  showing  rapid  needle 
lift  to  peak  pressure. 


frame  during  the  entire  film  sequence.  Timing 
marks  on  the  film  allow  estimates  of  the 
actual  film  speed  at  any  time. 

For  this  series  of  experiments  a  piston 
surface  simitar  In  general  shape  to  the  piston 
crowns  indicated  by  Kroeger  [  6  ]  and  by  Kato 
and  Onishi  [  7  ]  was  constructed.  In  order  to 
allow  optical  access,  a  truncated  piston 
surface,  as  shown  in  Fig.  4a  was  constructed. 
The  land  diameter  was  8.76  mm,  and  the  bowl 
diameter  was  87.3  mm  wfth  a  maximum  depth 
of  21.6  mm.  The  injector  was  placed  6.43  mm 
above  the  land  on  the  cusp,  providing  a  distance 
from  the  injsptor  tip  of  16  times  the  orifice 
diameter.  A  fiat  plate  was  attached  3.0  mm 
above  the  tipt  of  the  injector  to  simulate  me 
head  in  an  engine.  To  expose  the  spray,  the 
front  and  badk  of  the  piston  bowl  well  milled 
flat  to  the  Mjfiom  of  the  bowl,  leaving  a  center 
section  40.4  gmi  wide. 

The  spray;  waii  'mSk  lighted  using  me : 
copper  vapor  iaeer  beam  which  seta  expended 
and  coHlmaiai  through  a  set  of  leratas  and  then 
diffused  at  tita  entrance  window  of  me  bomb. 

m  iiw»  pnBib||n|jiw  ™  wwf  *  p»y» 

40. 


LON 

MED 

HIGH 

Rail  Pres. 
{MPa} 

4.7 

6.1 

7.8 

rMK  rT8S. 

(MPa) 

83 

91 

116 

Vol.  of  Fuel 
Inj.  (mm*) 

35 

69 

104 

Duration 

(ms) 

1.00 

1.92 

2.75 

Ave.VeL 

(m/s) 

267 

279 

292 

Table  1.  Experimental  Conditions 


The  resulting  films  were  analyzed  frame 
by  frame  using  a  512  by  512  CCD  video  camera 
connected  to  a  frame  grabber  board  in  a 
microcomputer.  The  digitized  frames  were 
processed  using  a  Sobei  edge  filter  [  10  ]  for 
finding  the  edge  of  me  spray,  a  2-D  median 
filter  to  reduce  image  noise,  and  thresholding 
to  improve  final  image  contrast.  The  result  of 
me  image  processing  is  a  two  level  image  of 
the  spray  and  chamber  edges.  The  spray  edges 
determined  by  this  procedure  were  then 
compared  to  the  parcel  location  pattern 
predictions  provided  by  the  model  at  the  same 
operating  conditions  for  the  same  time  from 
the  start  of  injection. 

RESULTS 

The  model  was  run  as  described  above  for 
the  set  of  conditions  used  in  the  constant 
volume  bomb.  A  typical  result  for  me 
'standard*  model  is  shown  in  Fig.  5a  where  me 
liquid  drop  parcel  locations  predicted  by  the 
model  are  indicated  by  the  circles.  Open  circles 
indicate  parcel*  whose  drops  have  not  hit  any 
•olid  surface.  The  parcel  size  Is  representative 
of  mo  rotative  size  of  the  drops  in  the  parcel, 
although  these  ore  rescaled  for  each  figure,  so 
size  evaluations  from  figure  to  figure  are  not 
JOtawye  valid.  The  SOM  Una  in  mo  figure 
indteitfes  the  outer  boundary  of  taaepray 
obtained  from  the  processed  high  speed  movies 
nt  the  conations  and  time  from  injection 
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Fig.  4a.  Top  and  aaction  views  of 
truncatod  piston  mad  In  photographic 
study. 


as  thft  numerical  prediction*,  in  this  figure  and 
the  following  figures,  the  radial  outline  of  the 
piston  bowl  is  shown  as  the  narrow  outline 
from  the  model  results,  while  the  thick  line 
indicates  the  truncated  bowl  outline  as  seen  by 
the  camera  (see  Fig.  4a).  The  spray  outline  is 
shown  for  the  frame  obtained  0.8  ms  after  the 
start  of  injection  for  the  fuel  rail  pressure  of 
6.1  MPa.  As  is  evident  from  the  figure,  this 
"standard”  model  does  not  accurately  predict 
the  general  spray  shape.  The  predicted  spray 
shows  little  dispersion  of  the  drops  relative  to 
the  photographic  results,  and  considerably  over 
predicts  the  penetration  of  the  spray  cloud. 
There  are  several  possible  explanations  for 
this  disagreement.  Most  arise  from  the  fact 
that  in  previous  work,  near  field  drop  behavior 
was  not  a  major  interest,  and  therefore  was 
not  considered  in  evaluating  the  accuracy  of 
the  results.  For  this  application,  near  field 
effects  will  clearly  have  a  major  impact. 

In  order  to  improve  the  results  of  the 
standard  model,  three  modifications  were 
proposed.  These  modifications  were  not 
intended  to  be  conclusive,  but  were  intended  to 
indicate  the  parametric  variation  and 
sensitivity  of  the  current  model  to  small 
modifications  in  order  to  improve  near  field 
predictions.  Conclusive  development  of  these 
submodel  improvements  will  require  more 
detailed  experimental  data  than  is  currently 
available. 
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The  first  of  the  proposed  modifications 
wit  to  altar  toe  size  of  tho  drops  as  they 
emerge  from  the  orifice  at  the  exit  of  the  fuel 
injector.  As  had  been  noted  [  2  ]  for  far  field 
comparisons,  choosing  a  single  initial  drop 
size  as  the  orifice  diameter  gave  reasonably 
good  agreement  for  downstream  drop  size  and 
velocity  for  dense,  solid  cone  sprays.  Since 
near  field  conditions  do  not  provide  the  length 
required  for  the  combination  of  drop  breakup 
and  coalescence  to  oblitoratp  the  initial  spray 
conditions,  initial  conditions  are  likely  to  be 
much  more  important.  For  this  numerical 
experiment,  an  initial  drop  diameter  an 
arbitrary  1/5  the  diameter  of  the  exit  orifice 
(80  |im)  was  chosen.  The  results  for  this 
calculation  are  shown  in  Fig.  5b.  The  model 
predictions  are  compared  with  the 
experimental  spray  outline  for  a  fuel  rail 
pressure  of  6.1  MPa  at  0.8  ms  aftsr  the  start  of 
injection.  This  figure  of  the  drop  parcel 
locations  still  does  not  agree  well  with  the 
spray  shape  from  the  photograph,  however 
some  efforts  not  seen  in  the  standard  case 
become  evident  These  include  a  number  of 
drops  which  are  entrained  in  the  induced  gas 
flow.  Because  of  their  small  size  and  low 


momentum,  these  drops  do  not  hit  the  spray 
impingement  surface.  The  overall  spray 
pattern  still  considerably  under  predicts  the 
lateral  spread  of  the  spray  cloud,  but  (foes 
come  much  closer  to  predicting  the  spray 
penetration  accurately. 

A  second  numerical  experiment  to  vary 
the  initial  drop  size  description  was  tested 
using  an  initial  drop  size  distribution  as 
proposed  in  the  original  KIVA  model  (Eqn.  1  ) 
with  an  average  Sautar  mean  dfometor  (SMO)  of 
50  pm.  The  results  for  this  calculation  are 
shown  in  Fig,  Sc.  As  with  the  previous 
numerical  esperimant  aH  other  contfitfons 
were  held  constant  In  eider  to  determine  the 
effect  of  toe  initial  drop  size  description  on 
the  overaft  spray  profile.  The  model  results  are 
rnmnyifl  mMi  hi  MMfimntfll  r—ulti  it  ft 

fuel  rail  pressure' of  6.1  MPa,  0.8  ms  after  the 
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While  it  is  quits  likely  that  a  more 
detailed  description  of  the  Injection  drop  size 
is  required,  it  is  equally  dear  that  modest 
modifications  of  this  parameter  alone  will  not 
produce  good  agreement  with  the  photographic 
record. 

A  second  set  of  modifications  to  the 
standard  model  was  a  set  of  modifications 
describing  the  change  in  the  spray 
characteristics  due  to  the  near  field  impact  of 
the  spray  on  the  solid  surface.  These 
modifications  included  changing  the  drop  size 
due  to  the  collision  with  the  solid  wail,  and 
modification  of  the  velocity  vector  of  toe  drop 
parcels  after  collision  with  the  wall.  The  first 
of  these  numerical  experiments  was  based  on 
the  proposal  that  in  toe  near  field  with  high 
velocity  liquid,  impact  with  a  solid  wall  is 
likely  to  have  a  major  effect  on  toe  drop  size 
distribution.  Previous  work  for  far  field  wall 
interactions  described  the  resulting  drop 
trajectory,  but  assumed  toe  drop  size  remained 
unchanged  in  toe  collision  [  5  ].  For  toe 
conditions  considered  in  that  problem,  the 
assumption  seemed  justified  by  reasonably 
good  agreement  with  some  quantitative  and 
some  qualitative  data.  For  the  high  velocity 
near  field  liquid,  it  seems  probable  that  the 
drop  size  distribution  is  significantly  altered 
by  the  impact,  generating  many  more  small 
drops  than  would  exist  without  the  impart. 

This  affect  was  qualitatively  indicated  by  the 
results  of  Wachters  and  Western ng  [11] 
who  showed  large  water  drops  at  approach 
Weber  numbers  of  184  break  up  into  many 
smaller  drops  at  a  surface,  in  an  attempt  to 
simulate  the  effect  of  change  in  drop  size 
distribution  due  to  wall  impact,  a  numerical 
experiment  was  run  using  the  standard  model 
description  for  the  initial  drop  sizes,  and 
allowing  all  drops  which  hit  the  surface  to 
Immediately  break  info  drops  1/5th  the  size  of 
the  drops  before  impingement,  with  mass  and 
kinetic  energy  conserved.  As  with  the  previous 
models  described,  toe  trajectories  of  these 
drops  were  determined  by  toe  jet  analogy 
model  of  Reitz  and  Naber.  As  seen  in  Fig.  5d, 
this  modification  does  have  a  significant 
impact  on  the  drop  parcel  distribution,  but 
when  compared  with  the  photographic  results, 
it  under  predfots  toe  rate  of  lateral  spread, 
while  predicting  the  penetration  depth 
reasonably  well. 


Fig.  5.  Comparison  of  the  expe/imental 
•pray  outlina  (bold  lino)  for  tho  modium 
rail  pressure  caso  (6.2  MPa)  at  0.8  ma 
•ftor  tho  start  of  Injoetlon  vs.  tho 
prodictions  of 

a. )  tho  "standard"  modal 

b. )  an  in)octod  drop  dia.  1/Sth 

orlflco  dia. 

c. )  an  Injected  drop  dia.  distribution 

(rn«25pm ) 

d. )  an  Implngod  drop  dia.  roducod  by 

1  /St  h. 


Tho  final  modification  applied  to  tho 
standard  model  was  to  incorporate  a 
component  of  velocity  normal  to  the  surface  in 
the  drop  trajectory  for  each  drop  that  hits  the 
surface.  The  drop  impingement  model  of  Reitz 
and  Naber  allows  drops  which  hit  the  surface 
to  leave  tangent  to  the  surface  in  a  direction 
relative  to  the  plane  of  incidence  calculated 
using  Eqns.  (  5  )  and  (  6  ).  For  near  field  high 
velocity  sprays,  it  is  likely  that  some 
component  of  velocity  normal  to  the  surface  is 
generated  by  the  impingement.  Detailed 
information  on  the  effects  of  drop 
impingement  at  high  Weber  number  (We-  pV2d/ 
o  )  conditions  is  scarce.  Wachters  and 
Westerling  [  11  )  provided  some  data  and  high 
speed  photographs  for  single  water  drops 
impinging  on  a  heated  surface  for  incoming 
Weber  numbers  up  to  184.  The  Weber  number  in 
the  current  set  of  experiments  is  much  higher, 
greater  than  1000.  Although  the  data 
presented  in  Wachters  and  Westerling  is 
primarily  for  heated  surfaces,  at  the  high 
Weber  numbers  cited,  the  heated  surface 
should  have  little  impact  for  drops  which  wet 
the  surface.  The  data  of  Wachters  and 
Westerling  indicated  that  for  high  incoming 
Weber  number  drops,  a  single  large  drop 
Immediately  shatters  into  many  smaller  drops. 
Their  photographs  (Fig.  11  of  (  11  ))  indicated 
that  many  of  these  smaller  drops  rebound  from 
the  surfaoe  with  a  small  velocity  normal  to  the 
surfaoe.  A  tummary  plot  from  Wachters  and 
Westsrilng  (  Fig.  7  of  [  11  p  suggests  that  for 
high  incoming  drop  Weber  numbers,  relatively 
tow  rebound  Weber  numbers  are  produced, 
however,  this  plot  represents  single  drop 
impingement  and  rebound,  without  subsequent 


breakup.  Although  the  data  of  Wachters  and  in  the  current  geometry.  As  for  the  low  fuel 
Westerling  is  for  incoming  Weber  numbers  injection  pressure  case,  the  spray  penetration 

much  lower  than  those  of  the  current  work,  the  is  somewhat  over  predicted  at  early  times 

qualitative  description  of  the  drop  breakup  and  after  the  start  of  injection.  Figure  9 
rebound  Is  accepted.  It  is  not  dear  exactly  summarizes  the  results  from  the  photographs 

what  form  the  normal  reflected  velocity  should  and  model  predictions  by  comparing  the 

have.  For  this  exploratory  study,  the  normal  calculated  radial  penetration  of  the  spray 

component  of  drop  velocity  after  impingement  verses  time  (indicated  by  the  lines)  for  the 

was  randomly  chosen  to  be  between  0  to  34%  three  rail  pressures  studied  (4.7,  6.1,  and  7.75 

of  the  incoming  velocity.  The  tangential  MPa)  with  the  results  from  the  photographs 

velocity  orientation  was  determined  as  before  (indicated  by  the  symbols), 
from  Eqns.  (  5  )  and  (  6  ).  The  drop  size  was  It  should  be  noted  that  there  remains 

reduced  to  1/5th  of  the  impinging  drop  size.  some  ambiguity  in  the  results,  due  to 

The  number  of  drops  in  the  parcel  was  chosen  characteristics  of  the  experimental  results 

to  conserve  mass  and  the  tangential  velocity  and  the  numerical  results.  The  experimental 

magnitude  was  chosen  to  conserve  kinetic  results  consist  of  the  high  speed  movies  of  the 

energy.  The  results  of  these  modifications  are  back-lit  spray  pattern.  The  start  of  injection 

shown  in  Figs.  6,  7,  8,  and  9.  These  results  could  not  be  determined  exactly  because  of  the 

represent  the  best  agreement  of  the  model  limited  frame  rate  of  the  camera.  For 

modifications  with  the  photographic  data  for  comparison  with  the  models  predictions  it  was 

the  limited  range  of  model  modifications  used.  assumed  injection  started  one  frame 

Figure  6  shows  the  set  of  five  spray  ( 0.2  ms )  before  the  first  appearance  of  the 

contours  taken  from  the  photographs  for  the  spray.  Since  the  pattern  is  so  dense,  spray 

fuel  rail  pressure  of  4.7  MPa  and  peak  injection  density  variations  are  not  clearly  seen.  The 

pressure  of  83  MPa.  The  figures  represent  photographic  results  show  a  generally  very 

times  from  the  start  of  injection  of  0.2,  0.4,  dense  cloud  of  spray  which  spreads  across  the 

0.6,  0.8,  and  1.0  ms.  The  model  predictions  for  piston  bore  with  time.  Since  the  photographs 

the  same  injector  conditions  and  time  represent  a  line-of-sight  average  attenuation 

intervals  are  superimposed  on  the  spray  of  the  incident  light,  the  visualized  "cloud* 

outlines.  The  agreement  between  the  two  is  may  be  a  uniform  cloud  with  the  outline 
reasonably  good,  although  at  early  times  after  indicated,  or  more  likely,  a  cloud  with  a 

the  start  of  injection  the  model  appears  to  thinner  center  section  as  indicated  by  the 

over  predict  the  spray  penetration.  Figure  7  model,  and  a  thick  outer  edge  where  the  slow 

shows  the  five  spray  contours  for  a  fuel  rail  drop  velocities  allow  for  increased  lateral 

pressure  of  6.1  MPa,  a  peak  injection  pressure  spread  of  the  spray  pattern.  For  sufficiently 

of  91  MPa,  for  the  same  time  intervals  as  in  dense  spray  clouds,  these  two  situations  are 

the  previous  test.  Again,  the  data  are  not  distinguishable  with  the  experimental 

compared  with  the  numerical  model  using  the  technique  currently  employed, 

modifications  for  surface  drop  breakup  and  It  should  also  be  emphasized,  that  in 

normal  velocity.  The  comparisons  again  show  applying  each  of  the  suggested  modifications 
relatively  good  agreement,  with  somewhat  to  the  standard  model,  some  arbitrary  choices 

better  agreement  for  the  spray  penetration  at  were  made  as  to  relative  initial  drop  size,  or 

early  times  in  the  injection.  Figure  8  shows  drop  size  after  surface  impact,  or  normal  drop 

the  five  spray  contours  for  the  fuel  rail  velocity  after  surface  impact.  There  exists 

pressure  of  7.8  MPa  and  peak  injection  very  little  data  under  the  conditions  of 

pressure  of  116  MPa.  Using  the  same  time  interest  for  detailed  analysis  of  these  aspects 

intervals  as  in  the  previous  figures,  the  data  of  the  problem.  One  objective  here  was  to 

are  compared  with  the  numerical  model  for  evaluate  which  mechanisms  appear  to  be  most 

these  input  conditions.  The  comparison  of  the  likely  to  play  a  major  role  in  the  problem  using 

experimental  spray  outline  and  the  predicted  some  reasonable  but  arbitrary  assumptions, 
parcel  locations  shows  reasonably  good  Subsequent  detailed  experiments  will 

agreement  between  fite  model  and  the  hopefully  allow  detailed  modeling  of  those 

experiment  for  this  set  of  operating  conditions  mechanisms  which  appear  to  be  most 

significant. 
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Fig.  6.  Comparison  of  experimental 
spray  outline  vs.  the  modified 
models  prediction  for  the  low  rail 
pressure  case  (4.7  MPa). 
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Fig.  7.  Comparison  of  experimental 
spray  outline  vs.  the  modified 
models  prediction  for  the  medium 
rail  pressure  case  (  6.1  MPa). 
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Fig.  8.  Comparison  of  experimental 
spray  outline  vs.  the  modified 
models  prediction  for  the  high  rail 
pressure  case  (  7.8  MPa). 
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Fig.  9.  Comparisons  of  ths  sxporimsntal 
impingod  sprsy  ponotratlon  vs.  ths 
modifiod  modols  prsdictions  at  ths 
thrso  rail  pressures  as  a  function  of  time. 

SUMMARY 


NOMENCLATURE 

d  orifice  diameter 

g(r)  drop  probability  distribution 

i  distance  from  nozzle  tip 

P  pressure 

r  drop  radius 

r3  2  Sauter  mean  drop  radius 

V  velocity 

a  impingement  angle  relative  to 

surface  normal 

p  probability  parameter 

p  density 

o  surface  tension 

v  kinematic  viscosity 

y  azimuthal  surface  angle 

subscripts 

g  gas 

rel  relative  |gas  -  drop| 


A  series  of  experiments  were  performed 
in  a  cold  high  pressure  constant  volume  bomb 
using  a  single  hole  fuel  injector  epray 
impinging  onto  a  land  in  the  middle  of  a  piston 
bowl.  The  overall  spray  pattern  was  recorded 
as  a  function  of  time,  and  compared  with  the 
results  of  a  numerical  spray  model.  After  the 
model  was  modified  by  including  estimates  of 
the  surface  interaction  effects  on  drop  size 
and  normal  velocity,  agreement  between  the 
experiment  and  model  was  reasonably  good. 

The  indication  from  this  result  is  that  the 
current  "standard"  modal  does  not  do  a  good  job 
of  predicting  spray  impingement  on  a  solid 
surface  near  the  fuel  injector.  Minor 
modification  of  the  current  model  improves 
this  agreement  Detailed  modeling  to  aooount 
for  this  improvement  awaits  more  detailed 
data  on  drop  sixes  and  velocities  near  solid 
aurfaoas  In  the  near  held  of  the  InioctiiiJ  sorav. 
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